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We have considered a cosmological model with a cosmological constant of the form A = 3a ^j +P% 7 a i P = const. 
The cosmological constant is found to decrease as t~ 2 and the rate of particle creation is smaller than the Steady 
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State value. We have found that this behavior gives -r 21 = 10 120 where Api is the value of A at Planck time. 
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fsji ■ Solutions with = 3a in the radiation dominated era and f3 = 6a in the matter dominated era are equivalent to the 
FRW results. We have found an inflationary solution of the de-Sitter type with (3 ~ 3 — 3a. Some problems of the 
Standard Model may be resolved with the presence of the above cosmological constant in the Einstein's equation. 

l/~} ' Since observations suggest a contribution of the vacuum energy density in the range 0.40 < O < 0.76, one gets 
4 < /3 < 12. If a = the minimum age of the universe is found to be H~ x (H p is the present Hubble constant) with 

\Q '. /?=oo. 
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1. Introduction 



Very recent results [1] suggest that the age of the universe is about 11 — 14 billion years. How- 
ever results obtained from Standard Model of cosmology are not consistent with these findings. 
Numerous models have been presented to reconcile these contradictions with observation. Only 
few of these models have given a satisfactory solution. The presence of a cosmological constant 
(A) in a given cosmology prolongs the age of the universe. The inflationary paradigm requires the 
Universe to have a critical density but observations do not support this. One of the motivations 
for introducing A term is to reconcile the age parameter and the density parameter of the uni- 
verse with current observational data [15, and references therein]. In an attempt to resolve these 
problems with current observations, we suggest a variation law for the cosmological constant of 
the form A = 3ot%2 + P% , at, j3 — const. Observational data indicate that A ~ 10~~ 55 cm -2 while 
particle physics prediction for A is greater than this value by a factor of order 10 120 . This discrep- 
ancy is known as the cosmological constant problem. It is interesting that this decay law helps to 
resolve this problem. The entropy problem that exists in the Standard Model can be solved by 
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the proposed decay law. The success of the inflationary model in extending the Big Bang models 
strengthen the theoretical prejudice toward a flat universe. When we allow G to vary with time we 
retain the conventional energy conservation law, i.e., the variation of A is cancelled by the variation 
of G. Though the proper way is to look for a field theoretic model for the variation of G and A 
nevertheless our present approach could be considered as a limiting case of some covariant theory 
yet to be discovered. Brans-Dicke theory [14] is one type of these theories but this theory allows 
only a decreasing G at the present time. The variation of G would have a significant effect on the 
evolution of stellar objects. Thus one way to obtain a restriction on the way in which G vary is 
via its impact on the evolution of astrophysical objects. 

A very recent discovery suggests that the universe is flat [12]. Some workers have 
shown that the universe is accelerating. In section 6. we show that the universe must be accel- 
erating if A > 0. This may be due to the fact that if gravity is increasing then the universe has 
to increase its expansion rate to escape the future collapse. Or alternatively, the decayed vacuum 
energy is given as a kinetic energy to accelerate the expansion of the universe. In this paper we will 
consider a flat universe which preserves the spirit of the inflationary scenario. The present model 
could resolve many of Standard Model problems with observation and thus could become a viable 
candidate as an alternative model. 

2. The model 

In a Robertson Walker metric, the Einstein's field equations with variable cosmological and gravi- 
tational 'constants' and a perfect fluid yield [2] 

3 T2 + W = 8*G P + A, (1) 

2 f + f + j|- 8 ^ + A - < 2 > 

where p is the fluid energy density and p its pressure. The equation of the state is usually given by 

p = ( 7 - l)p , (3) 

where 7 is a constant. Elimination of R gives 

3(p + p)R=-(*lp + p+J^)R. (4) 

3. Particle creation 

3.1 Matter- dominated universe(MDXJ) 

For a pressure-less MDU (p = 0) and for a constant gravitational constant (G), eq.(4) reads 

d(pR 3 ) _ R 3 dk 

dt ~ 8^G~dt ' ^ ' 



In this paper we will consider a decay law of the form [3,4] 

A = 3a| + 4, (6) 

where a and (3 are dimension-less constants. We suggest from the linear relationship between A 
and the Ricci scalar in the Einstein field equation that A has a general variation which resembles 
this scalar. For flat space one gets the above variation. For a flat universe (k — 0), eqs.(2) and (6) 
yield 

(2 - (3)RR = (3a - 1)R 2 , (7) 

which can be solved to give 

m = [ A(3 ~yy -^ (3 - 3 ^ , (8) 

where A = constant. Using eq.(8), eq.(6) becomes 

From eqs.(l), (6) and (8) the energy density takes the form 

p{t) = aJ^-^ ' (10) 

The vacuum energy density (p v ) is given by 



,,s A _ (2 - P)(6a - P) 1 

Pv{) 8nG 87rG(3-3a-^) 2 t 2 • l ' 



The deceleration parameter (q) is defined as 



RR 1 — 3a . . . . 

The density parameter of the universe (Q m ) is given by 

where p c = |^ is the critical energy density of the universe and H — j| is the Hubble constant. 
We notice that the Standard Model formula Vt m = 2q is now replaced by Q m = |g + |. However, 
both models give q = \ for a critical density. The density parameter due to vacuum contribution 
is defined as f2 A = -J^i- Using eqs.(8) and (9) this yields 

^g^f./*/'. (") 

We shall define f2 to tai as 

^total = fi™ + fi A , and ptotal = P + Pv ■ (15) 



Hence eqs.(13), (14) and (15) give $l to tai = ■'-■ This situation is favored by the inflationary scenario. 

_ (2-/3) x nm _2 (3-3a-/3) (6a - g) 2 

* p -(3-3a-/?)^ ' ^"3 (2-/3) ' ^ " "(2^-^ ' ^ 2 ' ^ 

(hereafter the subscript 'p' denotes the present value of the quantity). For ages larger than the 
Standard Model one requires (3 < 6a and for t p > 0, (3 < 2 and a > 1/3. This constraint indicates 
that A is positive. The precise value of a and (3 has to be determined from observational data. 
We now turn to calculate the rate of particle creation (annihilation) n, which is defined as [5] 

„=J-^!)| (17) 

Using eqs.(5), (8), (9) and (16) one obtains 

We remark that this rate is less than that of the Steady State model (= 3poH p ). li (3 = 6a then 
A = , n p = , t p = f-HT 1 and f2™ = 1. This case is equivalent to the Standard Model result. 
We observe that when a — |, q — and f3 decouples from all cosmo logical parameters. 

3.2 Radiation- dominated universe (RDU) 

This is characterized by the equation of the state p — |p (7 = 4/3). In this case, for a flat 
universe (k = 0), eqs.(l), (2) and (6) yield 

3(l-2a)|J + (3-2/3)| = 0. (19) 

This can be solved to give 

R = [ $-W) L U ](3-2/3)/2(3-3a-/3) £ = congt (20) 

(3 — 3a — p) 

For an expanding universe /3 < 3/2 and a > 1/2. Eqs.(l), (4), (6) and (20) give 

3(3-2/3)(3a-/3) l 3(3a-/3) 2 
4(3-3a-/3) 2 £ 2 (3-2/3) l J 

p = ^ ( 3 - 2 ^) 1 (22) 

^ 32vrG (3 - 3a - (3) t 2 ' V ' 

Thus an expanding universe ((3 < 3/2, a > 1/2) demands that the cosmological constant to be pos- 
itive. It is evident that if (3 — 3a the FRW results will be recovered viz. R oc i 1 / 2 , p oc t~ 2 , A = . 

4. No particle creation 

We now consider a model in which both G and A vary with time in such a way the usual en- 
ergy conservation law holds. 



4-1 Matter- dominated universe 

Equation (4) can be split to give [2,6] 

p + 3Hp = , (23) 

and 

A + 8nGp = . (24) 

Using eqs.(8) and (9), eqs.(23) and (24) yield 

R(t) = [ A Q-to-P) t] (2- m 3-3a-» > (25) 

where A = constant. 

(2-/?)(6a-/?) l 

A(t) = (3 -3a-/?)' ¥ ■ (26) 

p{t) = Fr 3(2-/?)/(3-3«-fl , F = const. (27) 

and 

U[) ~ [ AnF(3-3a-/3Y ' [ j 

Equation (24) represents a coupling between vacuum and gravity and that the vacuum decays 
to strengthen the gravitation interaction that will induce an acceleration of the expansion of the 
universe. Hence as long as gravity is increasing the expansion of the universe will continue. The 
variation of G could have been overwhelming in the early universe. This big gravitational force 
might have been the cause for stopping the rapid expansion during inflationary period and later 
assist in making the universe matter dominated. This because the increasing gravity forces smaller 
particles to form bigger ones. 

For (3 — 0, a — 0, G=const. and p = Dt~ 2 and R = [|At] 2//3 , which is the familiar FRW result. 
Moreover, the case (3 = 6a is equivalent to the Standard Model result. Clearly for j3 < 2 , a > 1/3 
the gravitational constant increases with time. In an earlier work [7] we have considered the effect 
of bulk viscosity in variable G and A models. We have shown that many of non-viscous models 
are equivalent to viscous models. The present model is equivalent to a viscous model with bulk 

viscosity (rj) varying as rj oc £ -n /( 1-n ) where r\ oc p n with n = 3 ( " 2 , . It has been shown that 
the development of the large-scale anisotropy is given by the ratio of the shear a to the volume 
expansion (6 = 3^) which evolves as [8] 

- oc £-( 3+3q - 2 ' 3 V( 3 - 3q -' 3 ) _ (29) 

9 

The present observed isotropy of the Universe requires this anisotropy to be decreasing as the uni- 
verse expands. Thus an increasing G guarantees an isotropic universe. For example, if 2/5 = 3 + 3a 
then G oc t -1 , R oc t 1 / 3 , p oc t _1 . This variation of G was considered by Dirac in his Large Number 
Hypothesis (LNH) model [9]. According to Dirac model one has a constant anisotropy and this 
may be a problem if the universe was not born isotropic because the anisotropy would not decay 
with time. 



4-2 Radiation- dominated universe 

Equation (4) now reads 

p + AHp = , (30) 

and 

A + 87rGp = 0. (31) 

Employing eqs.(20) and (21), eqs.(30) and (31) yield 

p = Ct -2(3-2/3)/(3-3a-/3) Q = [ ^ ~_ ^ _ U2(3«-ft/(3-3«-/9) ^ ^ = const . (32 ) 

We observe that when f3 = 3a the familiar FRW model is recovered. Abdel Rahman [2] has recently 
considered a closed universe model with a critical energy density where both G and A are variable. 
He found that R oc t ,Goct 2 , p oc t~ 4 in the radiation era. His solution corresponds to a — 1/2 
and a free (3. Thus both model, albeit different, evolve in the a similar way in the early universe. 
The large-scale anisotropy (see eq.(29)) becomes 

- OC £-3(3-2/3)/(3-3 a -/3) ^ p 3/2 (33) 

9 
Once again this anisotropy decreases with time as long as p decreases with expansion. 

4.3 Static solutions 

A static solution can be obtained for both matter and radiation dominated universes with (3 = 2 
and P = 3/2, respectively. Thus 

R = const., A = , p to tai = , n = . (34) 

It has been claimed by Kalligas et al. [10] that they have obtained a static universe with variable 
G and A. In fact, their solution is nothing but the above solution, since with R = const, eqs.(l) 
and (2) give A = so that G=const. Thus their claim of a static solution with variable G and A 
can not be true with p ^ —p. 

5. An inflationary solution 

This solution is obtained if we set H = const. Thus eqs.(l) and (2) give /3 = 3 — 3a so that 
A = 3H 2 . This can be integrated to give 



R = const exp(yA/3t) , ptotai — Pv , G = const. . (35) 

This is the familiar de-Sitter inflationary solution. 

6. An accelerating universe 

Now consider the case a — 0, i.e., A = /3^. Moreover, from eqs.(l)-(3) one finds 

A = 8j f (1 - b )pR + i R - (36) 



Thus writing A in the above form is equivalent to considering the universe to be filled with a fluid 
whose equation of state is given by p = — |p. Moreover, one may attribute the cosmic acceleration 
of the universe to the decay of the vacuum energy density (or A). It follows from eq.(8)-(14) that 
in the matter dominated Universe with G constant 

A((3-3) (0 _ 2]/(0 _ 3) /?(/?- 2) 1 J_(^-2)l_ 

K - [ (f3-2) tl ' A "(/?-3) 2 t 2 ' P ~ 47rG(/3-3)t 2 ' P * 6 > P * *' [6{) 

Pv ~ MG ~ MGj^WT^ ^ 3 - (38) 

* = -^ = ^ *> = g^fc 1 - "r = ig^l ^ 3 ' ^ 2 (39) 

^=(^2)^' n ? = 3(^' ^JTT)^ f* 2 (40) 

If G is changing with time one obtains from eqs.(28) and (29) 

G(t) = [ J^ZJL ]^/^), I oc i(3-^)/(^-3), /?^3, /?^2 (41) 

We observe that for Q^ > 0, /? > 2 and for fi™ > 0, (3 > 3. This implies that q p < and thus the 
universe must be accelerating at the present era. We also notice that G is an increasing function of 
time. Thus one may attribute the accelerating universe to the increasing gravity, i.e., the universe 
has to accelerate in order to overcome the ever increasing gravity. Others ascribed this to a scalar 
field that is still rolling down its potential "quintessence" to justify the acceleration of the universe 
[13]. We see that all the cosmological parameters shown above are functions of f3. The minimum 
age of the universe is H~ l and corresponds to f3 — oo. We restore the Einstein-de Sitter model 

when (3 — 0. We also note that ^- = (j§^u) 2 = 10 120 , where A Pi is the value of A at Planck 
time. This is evident from eqs.(9) and (21). This may help in solving the cosmological constant 
problem. Thus the cosmological constant is very small today because the universe is very old. We 
then have a model with one parameter that can be obtained from observation. Recent constrains 
on cosmological parameters from MAXIMA-1 [11] suggest that 0.40 < Q A < 0.76, this implies 
that 4 < j3 < 12. These results are found to be consistent with a flat universe. This constraint 
on (3 imposes a stringent constraint on all cosmological parameters so far known. The case (3 = 3 
corresponds to an inflationary solution in the matter dominated universe. If G is increasing with 
time the presently observed isotropy of the universe could be one of its consequences, as is evident 
from eq.(41). The case (3 = 2 represents a static universe, which is nonsensical for the present 
universe. 

7. Conclusion 

We have considered the decay law for A of the form A = Sa^ + (3^ , a, (3 = const. The cosmo- 
logical consequences of the model are shown to be very interesting. Many models in the literature 
can be retrieved from this model with particular choice of a and (3. The exact values of these 
constants still require more observational data to be determined. We note that the cosmological 



problems of the standard model, namely, the age problem, cosmological constant problem and the 
isotropy and homogeneity of universe, which are still unresolved could be tackled with our present 
approach. We have found that when a = the minimum age of the universe is H T 1 (corresponds 
to (5 — oo). As present estimates give Vt^ = 2/3 and Q™ = 1/3, we suggest a model of a = and 
(3 = 4, which gives an age of 2i/~ 1 and q p = —0.5, as describing our present universe. The infla- 
tionary solution, which solves the Standard Model problems, is shown to be built in the model. 
A model with a = gives rise to an accelerated expansion of the universe. This model is free 
from a lot of cosmological problems and could fit well with the present observational data. This 
model contains one free parameters viz., /3 that can be obtained from the present observational data. 
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